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Abstract 

We revisit the constraint on the primordial helium mass fraction Y p from obser- 
vations of cosmic microwave background (CMB) alone. By minimizing \ 2 of recent 
CMB experiments over 6 other cosmological parameters, we obtained rather weak 
constraints as 0.17 < Y p < 0.52 at la C.L. for a particular data set. We also study 
the future constraint on cosmological parameters when we take account of the pre- 
diction of the standard big bang nucleosynthesis (BBN) theory as a prior on the 
helium mass fraction where Y p can be fixed for a given energy density of baryon. 
We discuss the implications of the prediction of the standard BBN on the analysis 
of CMB. 



1 Introduction 



Recent precise cosmological observations such as WMAP pQ push us toward the era of 
so-called precision cosmology. In particular, the combination of the data from cosmic 
microwave background (CMB), large scale structure, type la supernovae and so on can 
severely constrain the cosmological parameters such as the energy density of baryon, cold 
dark matter and dark energy, the equation of state for dark energy, the Hubble parameter, 
the amplitude and the scale dependence of primordial fluctuation. 

Among the various cosmological parameters, the primordial helium mass fraction Y p is 
the one which has been mainly discussed in the context of big bang nucleosynthesis (BBN) 
but not that of CMB so far. One of the reason is that the primordial helium abundance has 
not been considered to be well constrained by observations of CMB since its effects on the 
CMB power spectrum is expected to be too small to be measured. However, since now we 
have very precise measurements of CMB, we may have a chance to constrain the primordial 
helium mass fraction from CMB observations. Since the primordial helium mass fraction 
can affect the number density of free electron in the course of the recombination history, 
the effects of Y p can be imprinted on the power spectrum of CMB. Recently, some works 
along this line have been done by two different groups El , which have discussed the 
constraints on Y p from current observations of CMB. In fact they claim different bounds 
on the primordial helium mass fraction, especially in terms of its uncertainty: the author 
of Ref. [2] obtained 0.160 < Y p < 0.501, on the other hand the authors of Ref. [3| got 
Y p = 0.250^q;q}4 at la confidence level. It should be noticed that the latter bound is much 
more severe than that of the former. If the helium mass fraction is severely constrained by 
CMB data, it means that the CMB power spectrum is sensitive to the values of Y p . In such 
a case, the prior on Y p should be important to constrain other cosmological parameters 
too and the usual fixing of Y p = 0.24 in CMB power spectrum calculations might not 
be a good assumption. Especially, analyses like Refs. jU 13 El predict light element 
abundances including 4 He from the baryon density which is obtained from the CMB data 
sets with the analysis fixing the value of Y p . Such procedure is only valid when Y p is not 
severely constrained by CMB. Thus it is very important to check the CMB bound on Y p . 

One of the main purpose of the present paper is that we revisit the constraint on 
Y p from observations of CMB alone with a different analysis method from Markov chain 
Monte Carlo (MCMC) technique which is widely used for the determination of cosmo- 
logical parameters and adopted in Refs. E] • in this paper, we calculate x 2 minimum 
as a function of Y p and derive constraints on Y p . We adopt the Brent method of the 
successive parabolic interpolation to minimize \ 2 varying 6 other cosmological parameters 
of the ACDM model with the power-law adiabatic primordial fluctuation. We obtain the 
constraint on Y p by this method and compare it with previously obtained results. 

We also study the constraint on Y p from future CMB experiment. A particular empha- 
sis is placed on investigating the role of the standard BBN theory. Since the primordial 
helium is synthesized in BBN, once the baryon-to-photon ratio is given, the value of Y p is 
fixed theoretically. Thus, using this relation between the baryon density and helium abun- 
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dance, we do not have to regard Y p as an independent free parameter when we analyze 
CMB data. We study how the standard BBN assumption on Y p affects the determination 
of other cosmological parameters in the future Planck experiment using the Fisher matrix 
analysis. 

The structure of this paper is as follows. In the next section, we briefly discuss the 
effects of the helium mass fraction on the CMB power spectrum, in particular its effects 
on the change of the structure of the acoustic peaks. Then we study the constraint on 
the primordial helium mass fraction from current observations of CMB using the data 
from WMAP, CBI, ACBAR and BOOMERANG. In section 4, we discuss the expected 
constraint on Y p from future CMB observation of Planck and also study how the standard 
BBN assumption on Y p can affect the constraints on cosmological parameters. The final 
section is devoted to the summary of this paper. 



In this section, we briefly discuss the effects of the change of the helium abundance on the 
CMB power spectrum. More detailed description of this issue can be found in Ref. |2j. 

The main effect of Y p on the CMB power spectrum comes from the diffusion damping 
at small scales. When Y p is large, since it is easier for electrons to recombine with 4 He than 
with H, the number of free electron becomes small. Thus the Compton mean free path 
becomes larger for larger Y p , which means that the diffusion length of photon becomes 
also larger. Since the photon-baryon tight coupling breaks down at the photon diffusion 
scales, the fluctuation of photon is exponentially damped due to the diffusive mixing and 
rescattering. Hence the CMB power spectrum is more damped for larger values of Y p . To 
see this tendency, we plot the CMB power spectra for several values of Y p in Fig. d We 
clearly see that the damping at the small scales is more significant for the cases with larger 
values of Y p . The effect of diffusion damping causes the change in the power spectrum at 
a percent level for 10 % change of Y p j^j. 

To see this more quantitatively, we consider the ratio of the second peak height to the 
first which is defined as JH] 



where (AT(/)) 2 = 1(1 + l)Ci/2ir. We do not discuss the first peak position and height 
because they are almost unaffected by the change in Y p . We calculate the responses of 
these quantities with respect to the change in the cosmological parameters around the 
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Figure 1: The CMB power spectra for the cases with Y p = 0.1 (green dashed line), 0.24 
(red solid line) and 0.5 (blue dotted line). Other cosmological parameters are taken to be 
the WMAP mean values for the power-law ACDM model. 



fiducial values, tu m = 0.U,u b = 0.024, fi A = 0.73, r = 0.166, n s = 0.99 and Y p = 0.24 
where o>j = Qih 2 with fi^ being the energy density of component i normalized by the 
critical energy density. The subscript b denotes baryon and m stands for matter which 
is the sum of baryon and CDM. h is the Hubble parameter, r is the reionization optical 
depth and n s is the scalar spectral index of primordial power spectrum. These values 
(except for Y p ) are the mean values of WMAP for the power-law ACDM model ^Hj. Also, 
we keep the flatness of the universe when we change parameters. We found 
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where we neglected the dependence on Q\ and r since their coefficients are very tiny even 
if compared to that of Y p . From these expressions, we can see that the response of C\ 
to the change in Y p is very sluggish. This is one of the reasons why we do not expect to 
obtain a meaningful constraint on Y p from CMB until recently. Moreover, the change of 
C\ caused by varying Y p is readily canceled by shifting other parameters. However, since 
observations of CMB now have become precise and cover wider multipole range, we may 
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have a chance to constrain Y p from current observations of CMB, which will be discussed 
in the next section. 

3 Constraint on Y p from current CMB observations 

Now we study the constraint on the helium abundance from current observations of CMB. 
For this purpose, we use the data from WMAP [I], CBI [TT] and ACBAR [12]. We also 
include the recent data from BOOMERANG experiment fSl El CHj To calculate y 2 
from WMAP data, we used the code provided by WMAP 021 EH EE! ■ For CBI, ACBAR 
and BOOMERANG, we made use of modules in COSMOMC [H]. As mentioned in the 
introduction, two groups have reported different bound on Y p using CMB data alone, 
especially in terms of its uncertainties. One group has obtained 0.160 < Y p < 0.501 [2] at 
la C.L., on the other hand the authors of Ref. [3] give the bound as Y p — 0.2501q q^4. The 
authors of Refs. [2] and [3J use the CMB data which cover the similar multipole region as 
that of ours. For details of their analysis, we refer the reader to Refs. [2] and [3J. If the 
severe bound on Y p is obtained from current CMB data, it means that the CMB power 
spectrum is sensitive to the value of Y p and the prior on Y p would affect the constraints on 
other cosmological parameters. Thus it is important to check the bound independently. 

For the analysis in this paper, we adopted a x 2 minimization by nested grid search 
instead of Markov chain Monte Carlo (MCMC) method which was used in their analysis. 
In our analysis, we apply the Brent method [20] of the successive parabolic interpolation 
to find a minimum with respect to one specific parameter with other parameters at a 
given grid, then we iteratively repeat the procedure to find the global minimum. For the 
detailed description of this method, we refer the readers to Ref. [21] . Here we assume a flat 
universe and the cosmological constant for dark energy. We also assume no contribution 
from gravity wave. In Fig. [2] we show the values of \ 2 minimum as a function of Y p . 
As seen from the figure, we do not have a severe constraint from current observations of 
CMB, which supports the result of Ref. [2]. Reading Y p values which give A% 2 = 1, we 
obtained the constraint at lcr C.L. as 0.17 < Y p < 0.52 for the case where the data from 
WMAP, CBI and ACBAR are used. When the data from BOOMERANG is added, we got 
0.25 < Y p < 0.54. We have also made the analysis for different data sets for comparison. 
In fact, we cannot obtain a significant constraint in the region 0.1 < Y p < 0.6 using 
WMAP data alone. Even if we add the data from BOOMERANG, we cannot constrain 
the value of Y p . Thus the data from CBI and ACBAR which cover high multipole regions 
are important to constrain Y p although the constraint is rather weak. 

As discussed in the previous section, the CMB power spectrum can be affected by 
changing the value of Y p . However this change can be canceled by tuning other cosmological 
parameters to give almost the same CMB power spectra. To see this clearly, in Fig. [3J we 
show the CMB power spectra for several values of Y p with other cosmological parameters 
being chosen to give almost indistinguishable angular power spectra. As seen from the 
figure, even if we take much larger or smaller values of Y p than usually assumed, we can 
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Figure 2: The values of A% 2 are shown as a function of Y p . Other cosmological parameters 
are taken to minimize x 2 . The red solid line is for the data of WMAP, CBI and ACBAR. 
The green dashed line includes BOOMERANG data in addition. 



fit such values of Y p to the data by tuning other cosmological parameters. 

When we include the data from BOOMERANG, the favored values of Y p are shifted 
to larger Y p . Notice that higher multipoles are more suppressed by increasing Y p , which 
is almost the same effect as decreasing rig^ 1 . Since BOOMERANG data favors red-tilted 
initial power spectrum compared to other data such as WMAP [22], it is reasonable that 
larger values of Y p are favored by BOOMERANG. Particularly, Y p = 0.24 which is used 
in usual analysis is just out of the la bound. However, it is allowed at 2a C.L. so we do 
not take this as a serious discrepancy from the standard assumption. 

4 Constraint on Y p from future CMB observations 
and the role of standard BBN theory 

In this section, we discuss the future constraint on the primordial helium mass fraction 
and other cosmological parameters. We especially want to investigate how the constraints 
are modified when we take account of the relation between ujb and Y p fixed by the standard 

# 1 Eqs. H2.3fl an d l|2.4|l show these properties quantitatively. We remark that they are also useful to 
understand the tendency that the values of n s which give the minimum x 2 f° r fixed Y p become larger 
as we increase Y p . This is because the suppression of higher multipoles caused by increasing Y p can be 
compensated by increasing n s . 
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Figure 3: The CMB power spectra for the case with Y p = 0.24, oj m = 0.135, Ub = 0.023, h = 
0.72, r = 0.117, n s = 0.96 (red solid line), Y p = 0.1,u m = 0.130, u b = 0.023, h = 0.72, r = 
0.101, n s = 0.95 (green dashed line) and Y p = 0.5, u m = 0.148, u b = 0.023, h = 0.72, r = 
0.117, n s = 0.99 (blue dotted line). Notice that these power spectra are almost indistin- 
guishable up to multipole region I ~ 1000. 



BBN theory. As mentioned in the introduction, the CMB anisotropies can be measured 
more precisely in the future, thus the primordial helium mass fraction may well be de- 
termined from CMB observations alone. The future constraints on Y p has already been 
investigated in Ref. j2] using the expected WMAP 4 year data, Planck and cosmic variance 
limited experiments. Here we study this issue supplementing the consideration regarding 
the prediction of the BBN theory. 

When we only consider observations of CMB alone, the primordial helium mass fraction 
Y p can be viewed as one of free independent parameters. However, when we take account 
of the BBN theory, Y p is not an independent parameter any more but is related to the 
value of the baryon-to-photon ratio or the energy density of baryon. Below, we discuss how 
such relation derived from the BBN calculation affects the determination of cosmological 
parameters in the future Planck experiment. 

First we give the relation between Y p and the baryon-to-photon ratio r\ from the calcu- 
lation of the standard BBN. rj and the baryon density ujb are related as 10 10 ?7 = 273.49cj&. 
Some groups have reported the fitting formula for Y p as a function of 77 |H1 123 EH EHj • 
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Here we adopt the fitting formula given in Ref. [H] 



ior P 



£ a nX n_1 + ^ b n x n ~ x AN + cnx 11 - 1 {AN) 2 + dnx 11 ' 1 {ANf 

n=l n=l n=l n=l 

x exp ^2 e nX n J , (4.5) 

where x = log 10 (10 10 r]), the coefficients a n , b n , c n , d n and e n are given in Ref. jH] and AN 
represents the number of effective degrees of freedom of extra relativistic particle species. 
The standard BBN case is obtained with AiV = 0. According to Ref. jH], the accuracy 
of this formula is better than 0.05 % for the range of 5.48 x 10~ 10 < 77 < 7.12 x 10~ 10 
(0.02 < Mb < 0.026) which corresponds to the 3a range obtained from WMAP and —3 < 
AiV < 3. Since scenarios with AiV ^ have been discussed in the literature including 
the possibility of negative AiV such as dark radiation in brane world scenario, varying 
gravitational constant and so on* 2 , we consider two cases when we discuss the future 
constraints. The first one is the case of the standard BBN, in other words, we assume that 
the energy density of extra radiation component as a fixed parameter with AiV = 0. For 
the other case, we treat AiV as a usual cosmological parameter which we vary, namely we 
assume AiV ^ 0. In this case, we use Eq. ()4.5|) to obtain Y p for given u b and AiV. 

Now we discuss the expected constraint from future CMB observations. For this pur- 
pose, we adopt the Fisher matrix method. Thus first we briefly review the Fisher ma- 
trix analysis which is widely used in the literature to study the future constraints on 
cosmological parameters. Detailed descriptions of this analysis method can be found in 
Refs. [HOI ED E21 EH] • For the CMB data, the Fisher matrix can be written as 

i x,y OXi OXj 

where X, Y — TT,TE, EE, Xi represents a cosmological parameter and Cov(C / x , Cj) 
is the covariance matrix of the estimator of the corresponding CMB power spectrum 
which is given explicitly in Ref. [31] • The la uncertainty can be estimated as (F _1 )jj 
for a cosmological parameter x,. For the fiducial model, we assumed the cosmological 
parameters as A = 0.86, uj m = 0.14,0^ = 0.024, = 0.73, r = 0.166 and n s = 0.99. Here, 
A represents the amplitude of scalar perturbation whose normalization is taken to be same 
as that of WMAP team ^U|. The fiducial value for Y p is fixed using the BBN relation 
Eq. ()4.5|) for a given ui b and AiV unless otherwise stated. A flat universe is assumed and 
we do not consider the contribution from the tensor mode. For dark energy, we assumed 

# 2 The negative values of AiV can also arise in a scenario with low reheating temperature T ro f ~ O(MeV). 
However, in this kind of scenario, the neutrino distribution functions are deviated from the thermal ones 
so the primordial helium abundance is modified in a way that the fitting formula Eq. (|4.5|l does not apply 
|2fil 1271 129| . The most recent Y p calculation in the low reheating scenario including the effects of 
neutrino oscillations are given in Ref. 29 . 
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Figure 4: Expected contour of the la constraint from Planck experiment. We take Y p 
as an independent free parameter. Other cosmological parameters are marginalized over 
in this figure. The thin and nearly horizontal band is the theoretical BBN calculation of 
Fj, as a function of Ub with its width representing la error from reaction rates. Here we 
assumed AiV = 0. 



the cosmological constant. To forecast uncertainties, we use the future data from Planck 
)3*4*] whose expected instrumental specifications can be found in Ref. j33| . 

Now we show our results. First we present the case with Y p being treated as a free 
parameter, namely we do not consider the BBN relation. In Fig. we show the expected 
la contour in the luj, vs. Y p plane. Other cosmological parameters are marginalized. We 
also draw the band for Y p as a function of Ub from theoretical calculation of the standard 
BBN with la error due to the uncertainties in the reaction rates. The uncertainty is 
dominated by that of the neutron lifetime, which is very small. As is clearly seen from the 
figure, the constraint on Y p from Planck is not significant compared to the uncertainties 
of the standard BBN calculation. Hence as far as we take account of the standard BBN, 
the helium mass fraction can be fixed using the BBN relation of Eq. (J4.5|) even for the 
precision measurements of CMB such as Planck. 

Next we discuss how the theoretical BBN relation can affect the determinations of 
cosmological parameters. First we consider the case with AiV = 0. When we determine 
the value of Y p for a given Ub using Eq. ()4.5|) (to be more specific, when we calculate 
numerical derivatives with respect to u>b, we simultaneously varied Y p following Eq. 1)4.50 ) 
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Table 1: Expected la uncertainties from Planck experiment using the temperature fluc- 
tuation alone. See the text for the fiducial values used in the analysis. 
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Table 2: Expected la uncertainties from Planck experiment using both the temperature 
and polarization data. See the text for the fiducial values used in the analysis. 



we can expect that the uncertainties of other cosmological parameters are reduced to some 
extent. In Tables E] and EJ we show the uncertainties of cosmological parameters from the 
future Planck experiment using the information of TT spectrum alone and that including 
polarization spectrum, respectively. The first and second rows in the tables correspond 
to the case without and with the BBN relation. For these cases, we assumed AiV = 0. 
Furthermore, we also show the cases with AiV ^ in the third and forth row in the tables. 

Now we discuss the cases with AiV = 0. As seen from the tables, when we assume 
the BBN relation, the uncertainties become smaller by a factor of 0(1) compared to that 
for the case with Y p being an independent free parameter. The parameter which receives 
the benefit most is n s . This is consistent with the fact that this parameter is the most 
degenerate parameter with Y p . Meanwhile, we note that the value of Y p fixed by the BBN 
relation for uj^ = 0.024 is slightly different from Y p = 0.24 which is usually used in the 
literature. We also evaluated the uncertainties fixing the helium mass fraction as Y p = 0.24 
independent of uib and checked that the la errors are quite similar to those for the case 
with Y p being related to u b by Eq. ()4.5|) . Since the change of C\ with respect to that of Y p 
is very small, in other words the derivative of C\ with respect to Y p is very small, we can 
have almost the same result even if we use a slightly different value for Y p . Thus we can 
just fix the value of Y p instead of using Eq. (j4.5|) even for the future CMB experiments 
such as Planck. 

Here the discussion for the case with AiV ^ is in order. Here we treat AiV as one of 
the cosmological parameters which should be varied. Notice that the addition of an extra 
radiation component can affect the CMB power spectrum through the speed up of the 
Hubble expansion and the early ISW effect because of the change of the radiation-matter 
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equality epoch. In fact, some authors have discussed the future constraints on cosmological 
parameters paying attention to AN [33 EE1 but without considering the effect of Y p . Here 
we study this issue allowing Y p to vary and also investigate the implications of the BBN 
relation on the future constraints on them. For the purpose of the Fisher matrix analysis, 
we assumed AiV = as the fiducial value. As already discussed, Y p , Ub and AiV are related 
by Eq. ()4.5|) from the BBN calculation. In Fig. we show the expected la contour in the 
cj& vs. AiV plane from Planck experiment for the cases with and without assuming the 
BBN relation. Naturally, when we assume the BBN relation, the uncertainties become 
smaller because it reduces the number of independent parameters. Also notice that, since 
the BBN theory relates Ub with AiV, the contour shrinks to the direction of correlation 
between Ub and AiV. 

In the third and forth rows of Tables ^ and the uncertainties of cosmological parame- 
ters are also shown for the cases with and without the BBN relation being imposed among 
Y p , u) b and AiV respectively. As expected, the uncertainties for cosmological parameters 
for the case with the BBN relation are smaller, however the differences are not so large. 

Here we comment on the implications of the BBN relation on the constraints from 
current observations of CMB. We have also made the analysis adopting the BBN relation 
to fix the value of Y p for given ujb- The constraint on a>& for this case is almost unchanged 
compared to the case with Y p = 0.24 being fixed. This can be more or less expected 
from the result we have shown in the previous section. This again shows that current 
observations of CMB are not so sensitive to the values of Y p . Thus predicting the helium 
abundance by the BBN theory using the CMB value of u>b, namely the procedure adopted 
in Refs. (1100111111, is valid at least with the current quality of the CMB data. 

5 Summary 

We revisited the constraint on the primordial helium mass fraction Y p from current obser- 
vations of CMB. Some authors have already studied the constraint [21 El, however their 
results were different especially in terms of the uncertainty. One of the main purpose of 
the present paper is to study the constraints on Y p from current observations adopting a 
different analysis method. Instead of MCMC method which was adopted by the authors 
of Refs. |2 El to obtain the constraint, here we adopted a \ 2 minimization by a nested 
grid search. We did not obtain a severe constraint in agreement with Ref. [2|. Using the 
data from WMAP, CBI and ACBAR as well as recent BOOMERANG data, we get la 
constraint as 0.25 < Y p < 0.54 and 0.17 < Y p < 0.52 for the cases with and without 
the data from BOOMERANG. It might be interesting to note that usual assumption of 
Y p = 0.24 is not in the la error range of BOOMERANG combined analysis but it is not 
of high significance at this stage so we can safely assume Y p = 0.24 for current CMB data 
analysis. 

We also studied the future constraint from CMB on Y p taking account of the stan- 
dard BBN prediction as a prior on Y p . Although we cannot obtain a severe constraint 
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Figure 5: Expected contours of la constraints from Planck experiment for the cases with- 
out the BBN relation (solid red line) and with it (dashed green line). We marginalized 
over other cosmological parameters in this figure. Y p is determined by the BBN relation. 



at present, observations of CMB can be much more precise in the future. Thus we may 
have a chance to obtain a precise measurement of Y p from upcoming CMB experiments. 
On the other hand, since the primordial helium has been formed during the time of BBN, 
once the baryon-to-photon ratio is given, the value of Y p can be evaluated theoretically 
assuming the standard BBN. Thus, in this case, we do not have to assume Y p as an inde- 
pendent free parameter when we analyze CMB data. We studied how such BBN theory 
prior on Y p affects the determination of other cosmological parameters in the future Planck 
experiment. We evaluated the uncertainties for the case with Y p being an independent free 
parameter and Y p being fixed for a given Ub using the BBN relation. We showed that the 
BBN prior improves the constraints on other cosmological parameters by a factor of 0(1) 
and also it induces some correlations among the parameters which appear in the BBN 
relation. As shown in Fig. HJ as far as we consider the standard scenario of cosmology, 
the helium mass fraction can be fixed for CMB analysis even in the future experiments 
since we can expect that the constraint from Planck is much weaker than the uncertainty 
of the theoretical calculation of the standard BBN. However, it is worthwhile to do CMB 
analysis treating Yp as a free parameter and measure the helium mass fraction indepen- 
dently from the baryon density since it provides a consistency test for the standard BBN 
theory (of course, measurements of primordial light element abundances by astrophysical 
means provide further consistency tests). By checking the robustness of the consistency 
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from various observations, the golden age of precision cosmology can push us toward the 
accurate understanding of the universe. 
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